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Abatraet.  Application  of  a  diagnoatic  cystem  to  a  hdicopter  gearbox  U  presented.  The  diagnostic  system  is  a 
nonparametiic  pattern  dassiiier  that  uses  a  malti-vahted  influence  mairix  (MVIM)  as  its  diagnostic  model  and 
benefits  from  a  fast  learning  algorithm  that  enables  it  to  estimate  its  diagnostic  model  from  a  small  number 
of  measnrement-fanlt  data.  To  test  this  diagnostic  system,  vibration  measurements  were  collected  from  a 
hdicopter  gearbox  test  stand  during  accelerated  &tigue  tests  and  at  various  fault  instances.  The  diagnostic 
results  indicate  that  the  MVIM  system  can  accurately  detect  and  diagnose  various  gearbox  faults  so  long  as 
they  are  included  in  training. 


1.  INTRODUCTION 

Helicopter  drive  trains  are  significant  contributors  to  both 
maintenance  cost  and  fiight  safety  incidents.  Drive  trains 
comprise  almost  30%  of  maintenance  costs  and  16%  of  me¬ 
chanically  related  malfunctions  that  often  result  in  the  loss 
of  aircraft  (Chin,  1993).  As  such,  it  is  crucial  that  faults 
be  diagnosed  in-flight  so  as  to  prevent  loss  of  lives. 

fhult  diagnosis  of  helicopter  power  trains  is  based  primar¬ 
ily  on  vibration  monitoring.  As  such,  considerable  effort 
has  been  directed  toward  the  identification  of  features  of 
vibration  that  are  affected  by  spedfic  faults  (e.g.,  Pratt, 
1986),  and  the  development  of  signal  processing  techniques 
that  can  quantify  such  features.  The  main  problem  with 
this  ^proacb,  however,  is  that  due  to  the  complexity  of 
helicopter  gearboxes  and  the  interaction  between  thdr  var¬ 
ious  components,  the  individual  vibration  features  do  not 
provide  a  reliable  basis  for  diagnoris. 

As  an  alternative  to  single-feature  based  diagnosis,  &ult 
diagnosis  based  on  several  features  can  be  performed  using 
pattern  classification  (Pau,  1981).  Among  the  various  pat¬ 
tern  classifiers  used  for  diagnosis,  artificial  neural  nets  are 
the  most  notable  due  to  their  nonparametric  nature  (inde¬ 
pendence  of  the  probabiUstic  structure  of  the  system)  and 
their  ability  to  generate  complex  decision  regions.  How¬ 
ever,  neural  nets  generally  require  extensive  training  to 
develop  the  decision  regions  (diagnostic  model).  In  cases 
such  as  helicopter  power  trmns,  where  adequate  data  may 
not  be  available  for  training,  artificial  neural  nets  may  mis¬ 
diagnose  the  fault. 

In  this  paper  we  demonstrate  the  application  of  a  diagnos¬ 
tic  method  that  can  establish  the  fault  signatures  based  on 
a  small  mimber  of  measurement-fault  drtta.  This  method 
utilizes  a  multi-valued  influence  matrix  (MVIM)  as  its  di¬ 
agnostic  model  which  provides  indices  for  diagnosability  of 
the  qrstem  and  variability  of  the  fault  signatures  (Daiuu 
and  Chin,  1991).  These  indices  are  used  as  feedback  to 
inqyrove  huilt  signatures  through  adaptation  (Clun  and 
Danai.  1993). 


2.  THE  MVIM  METHOD 

Measurements  are  processed  in  the  MVIM  method  as  il¬ 
lustrated  in  Fig.  1.  They  are  usually  filtered  first  to  ob¬ 
tain  a  vector  of  processed  measurements  P,  then  they  are 
converted  to  binary  numbers  through  a  flagging  operation 
(i.e.,  abnormal  measurements  characterized  by  ‘1’  and  nor¬ 
mal  ones  by  ‘0’),  and  finally  they  are  analyzed  through 
the  diagnostic  model.  In  the  MVIM  method,  flagpng  is 
performed  by  a  Flagging  Unit  that  is  tuned  according  to 
measures  of  diagnosability  and  fault  signature  variability 
obtained  fiom  MVIM  so  as  to  improve  the  fault  signa¬ 
tures.  The  MVIM  method  is  explained  in  detail  in  (Danai 
and  Chin,  1991)  and  (Chin  and  Danai,  1992).  We  w^  only 
discuss  the  overall  concept  here  for  completeness. 
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Fig.  1:  Processing  of  measurements  in  the  MVIM 
method 

2.1.  Fault  Signature  Representation 

Fault  signatures  in  the  MVIM  method  are  represented  by 
the  n  tmit-length  columns  Vy  6  TV"  of  a  multi-valued  in¬ 
fluence  matrix  (MVIM)  A: 

A  =  [V,,  ..,  V,-,  ..,  V,.  ]  (1) 

where  m  denotes  the  number  of  measurements,  and  n  rep¬ 
resents  the  number  of  faults.  Based  upon  this  influence 
matrix,  the  faults  can  be  ranked  according  to  their  pos¬ 
sibility  of  occurrence  by  the  closeness  of  their  influence 
vectors  Vj  to  the  vector  of  flagged  measurements  Y  (see 
Fig.  2).  In  the  MVIM  method  the  vector  of  diagnostic 
certainty  measures,  X,  which  ranks  the  faults  according  to 


I 


tlieir  possibility  of  occuzreace  is  defined  as 

i  ih  ■;  in)^  =  COs{oi,  Oj,  «»}’''  (2) 

where  the  Xj  teprescat  the  individual  iiagnostie  eertaini]i 
mtanms,  and  the  Oj  denote  the  individual  angles  between 
the  influence  vectors  V j  and  the  flagged  measurement  vec¬ 
tor  Y  (see  Rg.  2). 


total  number  of  false  alarms  and  undetected  faults,  and 
estimates  the  uniqtieness  and  variability  of  the  estimated 
fault  signatures  from  the  current  values  of  the  influence 
matrix,  so  thsd  it  can  use  these  measures  as  feedback  in 
the  next  ad^tation  roimd.  Adaptation  stops  when  the 
total  number  of  false  alarms  and  undetected  faults  are 
minimized,  and  the  uniqueness  and  consistency  of  fault 
signatures  are  enhanced. 


Y« 


Fig.  2:  Schematic  of  diagnostic  reasoning  in  the  MVIM 
method,  iUustrated  in  three  dimensional  space 

2.2.  Ertimation  gf  A 

The  influence  vectors  in  Eq.  (1)  are  not  known  arpriori  and 
need  to  be  estimated.  In  the  MVIM  method,  the  error  in 
dijigtinaia  is  Used  ss  the  basis  to  estimate/upda*;  the  in¬ 
fluence  vectors.  For  this  purpose,  the  &ult  signatures  are 
updated  recursively  after  the  occurrence  of  each  fault  to 
mmimize  the  som  of  the  squared  diagnostic  error  associ¬ 
ated  with  that  fault  (Danai  and  Chin,  1991). 

2.3.  Fault  Signature  Evaluation 

One  of  the  tmique  features  of  the  MVIM  method  is  its 
ability  to  evaluate  quantitatively  the  uniqueness  and  vari¬ 
ability  of  fault  signatures,  so  that  these  quantitative  mea¬ 
sures  can  be  used  to  improve  the  flagging  operation.  In 
the  MVIM  method,  the  uniqueness  of  fault  signatures  is 
represented  by  the  closeness  of  pairs  of  influence  vectors, 
and  the  index  of  ddagnosability  is  defined  to  characterize 
the  closest  pair  of  fault  signatures. 

In  the  MVIM  method,  the  variability  of  fault  tignatures 
is  defined  by  their  variance.  For  this  purpose,  the  vari¬ 
ance  matrix  of  A  is  estimated  to  provide  a  measure  of  the 
variations  of  individtial  components  of  A.  Since  in  the 
MVIM  method  the  components  of  A  are  adjusted  recur- 
rive^,  the  variance  matrix  can  be  readily  estimated  during 
training  (Chin  and  Danai,  1993).  The  index  of  fault  signa¬ 
ture  variability  is  defined  as  the  largest  component  of  the 
variance  matrix,  representing  the  largest  variability  in  the 
components  of  A. 

2.4.  Flagging 

The  influence  matrix  A  is  estimated  based  on  the  values  of 
the  flagged  measurement  vector  Y.  Thus,  before  the  influ¬ 
ence  matrix  is  used  for  diagnostic  reasoning,  the  integrity 
of  the  flagging  operation  needs  to  be  ensured.  Flagging  in 
the  MVIM  metlmd  is  performed  by  a  Flagging  Unit  that 
is  tuned  to  improve  the  diagnosability  of  the  system  and 
reduce  the  variability  of  the  fault  signatures  (Chin  and 
Danai,  1993). 

The  Flaggng  Unit  uses  a  sample  set  of  measurement-fault 
vectors  to  tune  its  parameters  iteratively.  After  each  pass 
throu^  the  training  batch,  the  Flagging  Unit  counts  the 


The  Flaggng  Unit  processes  the  residuals  as  follows  (see 
Fig.  3).  The  residual  vector  P  €  is  first  passed  through 
Hard-Limiter  I  (consisting  of  a  vector  of  m  thresholds,  hi,  , 
t  =  1, . . .  ,m)  to  produce  a  binary  vector  Z  €  This 
vector  is  than  multiplied  sequentially  by  the  normalized 
columns  of  a  Quantization  Matrix  and  then  thresholded 
by  Hard-Limiter  11  to  produce  the  individual  components 
of  the  flagged  measurement  vector  Y  €  B~.  Training  of 
the  Flagging  Unit  comprises  of  adjusting  the  thresholds  of 
Hard-Limiters  I  and  II  and  the  Quantization  Matrix,  as 
explained  in  detail  in  (Chin  and  Danai,  1992). 


Fig.  3:  Schematic  of  the  Flagging  Unit 


3.  EXPERIMENTAL 

Vibration  data  were  collected  at  NASA  Lewis  Research 
Center.  Various  component  failures  in  an  0H-5SA  main 
rotor  transmission  were  produced  during  the  experiments 
(Lewidd  et  al.,  1992).  The  configuration  of  the  trans¬ 
mission  which  was  tested  in  the  NASA  500-hp  Helicopter 
TVansmissioa  Test  Stand  is  shown  in  Fig.  4.  The  vibration 
signals  were  recorded  from  eight  piezoelectric  accelerome¬ 
ters  (frequency  range  of  up  to  10  KHz)  using  an  FM  tape 
recorder.  The  signals  were  recorded  once  every  hour  for 
about  one  to  two  minutes  per  recording  (at  the  tape  speed 
of  30  in/sec,  providing  a  bandwidth  of  20  KHz). 

A  total  of  five  tests  were  performed,  where  each  test  was 
run  between  nine  to  fifteen  days  for  ^proximately  four 
to  ei^t  hours  a  day.  New  components  were  used  at  the 
start  of  each  test.  When  a  component  fault  was  detected 
during  a  test,  it  was  replaced  with  a  new  one  for  the  re¬ 
mainder  of  the  test.  Eleven  failures  occurred  during  these 
tests.  The  status  of  various  faults  during  the  five  tests  are 
shown  in  Table  1  where  no-fault  cases  are  denoted  by  zo- 
In  Test  #1,  spiral  bevel  pinion  failure  ( ■^i)  is  estimated  to 
have  been  present  on  days  5-9  with  sun  gear  &ilure  (zi) 
also  occurring  on  day  9.  No  &ilures  occurred  in  Test  #2, 
so  all  the  nine  days  for  this  test  are  marked  as  normal. 
In  Test  #3,  planet  bearing  failure  (zj)  was  established 
to  have  been  present  on  days  3-4  and  11-12,  with  hous¬ 
ing  crack  (zs)  and  mast  bearing  foilure  (ze)  occurring  on 
day  9  and  day  13,  respectively.  Test  #4  is  estimated  to 
have  contained  planet  bearing  failure  on  days  11-12  and 
sun  gear  failure  on  days  14-15.  In  Test  #5,  housing  crack 
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4.  SIGNAL  PROCESSING 

In  order  to  identify  the  effect  of  faults  on  the  vibration 
data,  the  vibration  signals  obtained  from  the  five  tests  were 
digitized  and  processed  by  a  commercially  available  signal 
analyzer  (Stewart  Hughes,  1987)  with  four  processing  mod¬ 
ules:  (1)  Statistical  Analysis  (STAT),  (2)  Baseband  Power 
Spectrum  Analysis  (BBPS),  (3)  Bearing  Analysis  (BRGA), 
and  (4)  Signal  Averaging  Analysis  {SGAVJ. 

For  analysis  purposes,  only  one  data  record  per  day  was 
used  for  each  test.  The  data  records  were  taken  at  the 
beginning  of  the  day  unless  a  fault  was  reported.  When  a 
fault  was  detected,  the  record  was  taken  right  before  the 
fault  incident  to  ensure  that  the  data  record  reflected  the 
fault.  Also,  in  order  to  reduce  estimation  errors,  each  data 
record  was  partitioned  into  sixteen  segments,  and  param- 
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Fig.  5;  Configuration  of  the  MYIM  system  as  apphed  to 
the  OH-58A  main  rotor  transmission. 

#3  and  #4  contained  most  of  the  ftulure  modes  (i.e.,  4 
out  of  6).  Therefore,  the  parameters  horn  these  two  tests 
were  used  to  train  the  MVIMs.  Not*  that  not  all  of  the 
failure  modes  were  included  in  training,  so  the  test  results 
were  not  expected  to  be  perfect.  For  training  the  detection 
MVIMs,  only  the  19  parameters  from  the  STAT,  BBPS, 
and  BRGA  modules  were  used.  Previous  studies  on  this 
data  show  that  these  19  parameters  are  adequate  for  de¬ 
tection  (Chin,  1993).  For  training  the  diagnostic  MVTMs, 
all  of  the  54  parameters  were  utilized. 

The  initial  values  of  the  detection  MVIMs  ( 19  x  2)  and  di¬ 
agnostic  MVIMs  (54  X  7)  were  each  set  to  0  (i.e.,  matrices 
with  all  zero  entries),  and  the  initial  values  of  the  Quan¬ 
tization  Matrices  (19  x  19  for  detection)  and  (54  x  54  for 
diagnosis)  were  set  to  identity  matrices.  The  initial  thresh¬ 
old  levels  for  Hard-Limiter  I  were  set  at  the  mean  plus  one 
standard  deviation  of  the  corresponding  parameter,  and 
for  Hard-Limiter  II,  they  were  set  at  0.5.  The  maximum 
number  of  epochs  for  training  the  detection  and  diagnos¬ 
tic  MVIMs  was  set  to  50.  After  each  epoch,  the  detec¬ 
tion/  diagnostic  performance  of  MVIMs  within  the  training 


set  was  tested,  l^aiaing  was  stopped  once  perfect  detec- 
tioa/diagnosis  was  achieved,  to  avoid  overtraining  (Hertz 
et  at,  1991).  The  number  of  epochs  used  for  individual 
detection  MVIMs  were:  8,  5,  SO,  37,  50,  15,  50,  and  50 
for  aocelerometem  #1  to  #8,  respectively,  whereas  for  di¬ 
agnostic  MVIMs  they  were:  50,  1,  2,  2,  26,  50,  50,  and 
SO.  According  to  the  number  of  epochs  used  for  individual 
MVIMs,  it  is  clear  that  the  detection  MVIMs  assodated 
with  accelerometers  #3,  #5,  #7,  and  #8  did  not  achieve 
perfect  detection  within  the  training  set.  Similarly,  the  di¬ 
agnostic  MVIMs  associated  with  accelerometers  #1,  #6, 
#7,  and  #8  did  not  achieve  perfect  diagnosis  within  the 
training  set. 

The  MVIMs  trained  on  Tests  #3  and  #4  were  evaluated 
for  all  of  the  tests.  For  this  purpose,  the  19  parameters 
horn  each  of  the  eight  accelerometers  were  first  passed 
through  the  corresponding  detection  MVIM  for  all  of  the 
five  tests  to  reflect  the  occurrence  of  faults.  Once  a  fault 
was  posted  by  a  detection  MVIM,  the  set  of  54  parame¬ 
ters  from  that  accelerometer  was  passed  through  the  cor¬ 
responding  diagnostic  MVIM.  At  the  final  stage,  the  diag¬ 
nostic  certainty  measures  obtained  from  the  eight  diagnos¬ 
tic  MVIMs  were  consolidated  by  the  voting  sdieme,  which 
utilized  weights  reflecting  the  speed  of  convergence  of  in¬ 
dividual  fault  signatures  during  training  (Chin,  1993).  As 
such  those  influence  vectors  which  converged  faster  were 
assigned  larger  weights  and  vice  versa.  Zero  waghts  were 
asagned  to  the  influence  vectors  which  did  not  converge 
during  training,  and  unity  weights  to  those  which  con¬ 
verged  within  one  epoch. 

The  diagnostic  results  obtained  from  the  above  diagnostic 
system  ate  shown  in  Table  2.  The  results  indicate  that 
the  MVIM  system  was  able  to  produce  perfect  diagnostics 
for  Tests  #3  and  #4,  on  which  it  was  trained,  and  that  it 
diagnosed  88%  of  the  faults  in  all  of  the  tests.  Specifically, 
the  results  in  Table  2  indicate  that  the  MVIM  system  pro¬ 
duced  two  false  alarms  (on  day  4  of  Test  #1  and  day  6  of 
Test  #5),  and  five  misdiagnoses  (on  days  5-8  of  Test  #1 
and  day  11  of  Test  #5).  In  addition,  this  system  produced 
equal  diagnostic  certainty  measures  for  the  no-fault  case 
and  sun  gear  failure  on  day  10  of  Test  #5  and  could  only 
diagnose  one  of  the  faults  on  day  9  of  Test  #1,  and  on 
days  10  and  11  of  Test  #5.  However,  it  should  be  noted 
that  faults  and  zj  were  not  included  in  training,  so 
no  fault  agnatuies  were  estimated  for  them.  The  correct 
diagnostic  rate  of  MVIM,  with  these  two  faults  excluded 
would  be  over  95%,  which  is  quite  noteworthy  conadering 
that  the  MVIM  system  was  trained  on  a  small  number  of 
measurement-fault  data  with  very  few  repetitions  of  each 
&ult. 

6.  CONCLUSION 

A  fault  diagnostic  system  based  on  the  MVIM  method  is 
iq>plied'to  a  helicopter  gearbox  for  which  only  a  small  num¬ 
ber  of  measurement-fsult  data  was  available.  This  diag¬ 
nostic  system  utilizes  two  levels  of  isolation  and  integrates 
the  results  obtained  from  various  accelerometers  throu^ 
a  voting  scheme.  The  diagnostic  results  indicate  that  the 
MVIM  system  correctly  detected  all  of  the  fault  incidents, 
with  only  two  frJse  alarms  posted.  The  results  further  in¬ 
dicate  that  this  system  correctly  diagnosed  all  of  the  faults 
it  was  trained  for. 


Table  i:  Estimated  faxJt  stains  from  each  day  of  ike  five 
tests.  Same  notations  are  adopted  as  >n  Table  I 
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